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Abstract 

Using the quantum magnetohydrodynamics (QMHD) model, linear dispersion of magnetosonic 
(— I ■ waves are studied in a quasineutral quantum electron-ion plasma in two distinct regimes of nonrel- 

ativistic and relativistic degeneracies considering also the plasma composition effect. It is shown 
that the paramagnetic spin effects of the degenerated electrons plays a key role in dynamics of 
"q^ magnetosonic waves. The linear wave-speed is found to have minimum value at some degeneracy 

parameter in such plasmas. This is due to delicate interplay between relativistic degeneracy and 
the Pauli spin-magnetization. It is also revealed that the plasma composition has significant effect 
Q_J on the linear dispersion in the relativistic degeneracy limit contrary to that of nonrelativistic case 

and Zeeman energy has significant effect in nonrelativistic degeneracy regime unlike that of rela- 
tivistic one in the linear perturbation-limit. Current findings can have important applications in 
both inertial plasma confinement and astrophysical degenerate plasmas. 

PACS numbers: 52.30.Ex, 52.35.-g, 52.35.Fp, 52.35.Mw 
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I. INTRODUCTION 



A quantum plasma is usually referred to as a degenerate dense and cold ionized matter. 
Despite the high temperatures observed for astrophysical compact objects, it has been shown 
l| that these objects can be regarded as fully degenerate cold quantum plasmas. This is 



because the temperature of plasma species in these objects stay below a well defined 
temperature of order lO'^K. Vast majority of the recent studies on quantum plasmas 



ermi 



m 



have been inspired by the pioneering works of Bohm, Pines and Levine et.al |12|-|14| and the 
earlier works of Chandrasekhar. Distinction between an ordinary and quantum plasma (also 
known as the zero-temperature Fermi- plasma) is that, in a quantum plasma contrary to the 
ordinary one, the inter-fermion distances are much lower than the characteristic de Broglie 
thermal wavelength, h/ {2TimkBTY^'^ so that the particle waves overlap causing a new type 
collective phenomenon through the electron degeneracy and tunneling effects. These prop- 
erties which are encountered only in quantum plasmas, for instance in ordinary metals and 
compact astrophysical objects, lead to distinct features of wave dynamics when the electron 
number-density or the ambient magnetic field strength or even the system dimensionality 
and the degrees of freedom of fermions is varied ISMlTj]. The electron degeneracy is a di- 
rect consequence of the Pauli exclusion principle and effectively rules the thermodynamical 



properties of a dense matter 



18|- 



Most of current investigation on quantum plasmas are motivated due to the applicability 
of the quantum magnetohydrodynamics (QMHD) to semiconductors [191]. Quantum opti- 
cal and electronic transport phenomena are also known as main ingredients in the study 
of metallic and semiconductor nano-structured materials such as, nano-particles, quantum- 



wells, quantum-wires and quantum-dots 



20l |. Therefore, the subject of dense cold ionized 



matter under strong magnetic field can promise wide applications in both 



aboratory as well 



as naturally occurring plasmas such as astrophysical compact objects [21|, |22|. One of the 
differences of the QMHD model from those of classical counterparts is the inclusion Bohm 
tunneling effect. It has been shown that the increase in the electron number density enhances 
;his feature significantly altering the dynamics of linear and nonlinear wave propagations 



231]. More recently, Marlkund and Brodin have extended QMHD model and incorporated 



the electron spin-1/2 effects 2J] by introducing a generalized term for the quantum force a 
feature which is expected to be dominant in superdense plasmas such as Fermi-Dirac. It has 



2 



been estimated that relativistic electrons in a white dwarf can generate magnetic fields of 



the order of 10^ G and higher 25|, and that in a neutron star the self- field may be of order 



10^'^ — lO^^G [26|. On the other hand, it has been shown that, in a perfectly conductive 
quantum plasma, the spin magnetization term introduced in quantum magnet ohydro dy- 
namics (QMHD) equations acts as a negative pressure-like entity significantly affecting the 



dynamics of spin- induced nonlinear magnetosonic excitations 27 1. 



In 1935 Chandrasekhar 



in a pioneering work, proved that the gigantic inward grav- 



itational pressure can lead to the new state of relativistic degeneracy for electrons and may 



cause the ultimate collapse of the star in extreme degeneracy cases [29| . The sudden collapse 
of the white dwarfs was shown to be due to softening of the degeneracy pressure of electrons 
relativistic motion. This effect has been shown to alter the whole thermodynamical proper- 
ties of the star [sol . Current investigation is an attempt to trace possible differences in linear 
magnetosonic wave dynamics caused due to change in the plasma degeneracy state such as 
nonrelativistic and relativistic ones and explore the composition effects on linear dispersion. 
The presentation of the article is as follows. The QMHD plasma model including the spin 
contribution is introduced in Sec. [Tll The linear dispersion of spin-induced magnetoacoustic 
waves are derived in Sec. Illll for different plasma degeneracy regimes and the numerical 
results are presented in Sec. [IVl 

II. QMHD MODEL INCLUDING SPIN EFFECTS 

Let us consider a collisionless dense degenerate quasineutral electron-ion perfectly con- 
ducting plasma with electrons and ions immersed in a spatially varying magnetic field in a 
fixed direction. We further assume that the magnitude of the external field is very high and 
the dynamics of wave propagation is strongly affected by spin contribution from fermions. 
In the center of mass frame of plasma the continuity relation reads as 

|^ + V-(pV) = 0, (1) 

where, p = merie+mirii ~ mjnj and V = {riemeVe+nimiVi) / p are the plasma center of mass 
density and velocity. The momentum equation, taking into account the quantum diffraction 



3|] and spin pressure [2J] effects, in the center of mass frame takes the form 27|, l31| 
dV 

_ + (V ■ V) V = ( J X B - VP + Fq) , (2) 



where, the quantum force Fq = Fb + Fs is a collective contribution of quantum Bohm-force, 
Fb, and the so-called spin-force, Fs, namely. 



'Q 



-V^^ + MVB. 



(3) 



2memj ^p. 

The quantities J, B and P are the magnetization current, magnetic field and scalar total 
center of mass plasma pressure, respectively. The magnetization current, J, is related to the 
spin-magnetization (Pauli-magnetization) per unit volume Mp = {?)^\ne/2kBTFe) B {Tpe 
being the electron Fermi-temperature) through; J = Hq^V x (B — poM), where, hb = 
eh/2meC is the Bohr magneton and h is the normalized Plank constant. It should be noted 
that, in the case of fully degenerate quantum plasma the Pauli paramagnetism dominates 
the Langevin-Type one, hence, the Fermi-temperature, Tpe, replaces the thermal electron 
temperature, Te in magnetic pressure term. In a completely degenerate plasma, due to the 
Pauli exclusion principal, only a small fraction (Tg/Tpe), i-e. unpaired free-electrons, with 
energies nearly equal to the Fermi-energy contribute to the collective spin effects. This effect 
gives rise to the well-known Pauli paramagnetic law of the form mp = ?)N / {2Epe) 18| 
instead of the classical law of nic = NfiB tanh.{ fiBB /kBTe) frequently used for dense plasmas 
in literature. On the other hand, the Landau demagnetization factor of 1/3 of this value 
should be subtracted, which results in total sum of M = Mp + = ne^%B/ Ep^. The 
susceptibility of weakly interactin g re lativistic Fermi-gas in a weak magnetic field [fipBo -C 
nieC^) has been evaluated in Ref. [3^, which is as follows 



Xp 





(4) 



where, 7 and Ppe are the s-wave interaction parameter and the Fermi relativistic momentum, 
respectively. It can be shown that in the nonrelativistic limit with 7 = the susceptibility 
given in Eq. (j3j) reduces to the well-known Pauli paramagnetic susceptibility. In order to 
complete the set of QMHD fluid equations, we use the generalized Faraday law without the 
Hall term in the form 



dB 

'dt 



V X <^ V X B - r^J 



J X B + Fq 


irte 


' d 


/V X B\ 






di 


V epone / 



■ V 



V X B 



(5) 



where, 77 is the plasma resistivity. Moreover, with p = p{x, t) and V = V{x, t), and assuming 
the propagation of the magnetoacoustic wave to be in x direction perpendicular to the 
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magnetic field B{x,t) along the z-axis, we may write simplified one-dimensional version of 
the complete QMHD equations as j^^] 



|f + ^ = 0, p = p{x,t), V = Vix,t), B = Bix,t), 



dV _|_ ydV _|_ _B_dB_ _|_ 1 aP _ _ji_d_ ( J_^VpA 

dt dx flop dx p dx 2memi dx \ ,Jp dx"^ j 



4 

meTnipc^ dx ^ \ Pp^ j \/ ^ ' V "leC 



ae I dBv _ rj_a'^B ^ g 

dt dx Ho dx'^ 



(6) 



On the other hand, in a relativistically degenerate zero-temperature Fermi-Dirac plasma one 
can ignore the ion-pressure compared to that of the dominant electron degeneracy pressure 
33 1 which can be expressed in the following general form [l| 



= |i? (2i?2 _ 3) v^rT^ + 31n [r + VlTW] } , (7) 

in which R = (Ppe) / {rriec) = {p/PcY^^ {pc — 1-97 x lO^gr/cm^) is the relativity parameter 
with Ppe being the electron Fermi relativistic momentum. The relativistic degeneracy due 
to the Pauli exclusion principle increases with decrease in inter-fermion distances in degen- 
erated plasmas. It is noted that, although the relativistic effects arise due to increase in 
fermion number- density in a dense degenerate plasma, however, unlike for the low-pressure 
relativistic plasmas the degeneracy pressure in the fermion momentum fluid equation usually 
dominates the relativistic (electron/ion) dynamic effects in super-dense plasma state and the 
relativistic electron/ion mass effects may be ignored. It has been shown 2^ that for a zero- 
temperature dense degenerate Fermi-gas such as a white-dwarf with a mass-density, p, the 
degeneracy pressure turns from oc p^^^ (with polytropic index 3) dependence for normal 
degeneracy in the limit i? ^ 1 to Pg oc p^/^ (with polytropic index 2/3) dependence for 
relativistic degeneracy case in the limit R^ 1. On the other hand, the relativistic degener- 
acy starts at mass density of about 4.19 x 10^{gr/cm^) of the order in the core of a 0.3Mq 
white dwarf, which corresponds to a Fermi- momentum Ppg ~ 1.29meC (equivalent to the 
relativistic degeneracy parameter value of i? ~ 1.29) or the threshold velocity of upe ~ 0.63c 
(the Fermi relativistic factor ~ 1.287). This may justify our nonrelativistic treatment 
of electron mass effects in above arguments. 



The 
et al. 



Droblem of spin-induced magnetosonic has been thoroughly investigated by Masood 



34l | for the case of the relativistic degeneracy, however, for completeness in cur- 



rent work we try to extend the previous study to the case of quantum plasma with the 



Coulomb interaction effect. In a dense quantum plasma the coupling parameter T repre- 
senting the potential to kinetic energy ratio may exceed unity leading to strong coupling due 
to other many-body effects such as Thomas-Fermi screening, Coulomb interactions, electron- 
exchange and ion-correlation. Salpeter 35| has examined these effects in the relativistically 
degenerate plasmas and concluded that the first-order correction to the degeneracy pressure 
can be that of Coulomb interactions and other effects may be negligible. Following the 
Salpeter expression for classical Coulomb plus Thomas-Fermi corrections, we may write 



C+TF 



^^2/3 / 4 \ 1/3 i^aZ^/^f f aV'"" 



107r2 VgvrJ ^^ 175 97r2 \9ti) yJl + R^ 



(8) 



where, a = e^/hc ~ 1/137 and Z are the fine-structure constant and the atomic number, 
respectively. The first term corresponds to the Coulomb negative pressure and the second 
term (proportional to a^) is a relatively smaller contribution due to Thomas- Fermi screening, 
which will be neglected in current scheme. 

On the other hand, the plasma resistivity, i], of a completely degenerate Fermi-Dirac 
plasma is supposed to be negligible, since, the electron-ion collisions are lowered by Pauli- 
blocking mechanism, i.e., ~ 0. By using Eq. ([6]) this lead to the simple relation, B/Bq = 
p/po = p (the so-called frozen-in field condition). The new quantities Bq, po and p denote 
the equilibrium values of magnetic field intensity, plasma mass-density and the normalized 
mass density, respectively. Therefore, the complete set of dimensionless equations governing 
the dynamics of magnetosonic waves reduce to 

dp I dpV _ ri 

at dx 

f + ff= (9) 



dx dx 



where, /3 = a{2^/^ /5){3Z^ /ny/^ and we have used (l/p)9^.Pe(p) = d^VTTW and 
(1/ p)dxPc{.p) = dx{/3R) and employed the following scahngs in order to obtain the di- 
mensionless equations ignoring the normalization bar notation 

X —X, t -)■ — , p ppo, y yVs, (10) 

where, c, Upi = ■>/ e^rieo/ {sonii) and Vs = c^jmjnii are the vacuum light speed, the char- 
acteristic ion plasma frequency and ion sound-speed (this speed is much higher despite the 
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name comparable to the Fermi-speed of an electron in a solid), respectively, and the parame- 
ter po denotes the equilibrium plasma mass-density. We have also introduced new fractional 
plasma entities used frequently in forthcoming discission, namely, the quantum diffraction 
)arameter, H = rrii /2me{hujpi) / {rrieC^) , relativistic degeneracy parameter, _Ro = (Po/PcY^^ 



331], and the normalized Zeeman energy, e = ^bBq/ {meC 
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III. COULOMB INSTABILITY AND DEGENERACY REGIMES 

Now let us consider the special cases of ^ 1 and ^ lin the spin-contributed term in 
the momentum equation (Eq. Q). These limits will be shown to lead to distinct features in 
dispersion of linear spin-induced magnetosonic waves. Since, the normalized Zeeman energy 
parameter is solely defined by the external magnetic field strength, the two extreme degen- 
eracy limits, namely, the nonrelativistic {Rq ^ 1) and relativistic {Rq ^ 1) degeneracy can 
be distinguished for a given magnetic field strength value, say, Bq = lO^T (e ~ 10~^) typical 
of spin-dominated compact astrophysical objects. On the other hand, the reader should 
have already realized the elegant relationship between the quantum diffraction parameter, 
H, which is due to the electron tunneling and the relativistic degeneracy parameter, Rq, due 
to the increase in the electron number density, namely, H = eh\/ pcR^jii ji^mj'^c?). Fourier 
analyzing Eq. ([9]) results in the following dispersion relation for the linear magnetosonic 
waves in the coupled relativistically degenerate plasma case 



+ 



A;2 2 3 



R. _^ 



+ (11) 



The differences in dispersion relation (Eq. f|TT]) ) between the two cases of -C 1 (relativistic 
degeneracy) and ^ 1 (normal degeneracy) is readily apparent from the dispersion relation 
Eq. ( ITT]) . It is observed that, for the case of relativistic degeneracy the last term in dispersion 
relation is negligibly small and the phase speed of the linear magnetosonic-wave is almost 
independent from Zeeman energy, e, while for the nonrelativistic plasma regime this is not 
the case. On the other hand, for the case of nonrelativistic degeneracy the first term in the 
dispersion relation is insignificant and the wave becomes dispersionless. The approximate 
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dispersion regimes in extreme degeneracy limits can be summarized as below 



2 _ IP^ , 



fc2 ~ ^2 -r 3 



Ro 



-/3 



< 1 



-/3 



(12) 



It is apparent from Eq. (fTTl) that the instability can set in when 



H^e Ro 

h — 

2 3 



-(3 



(13) 



This condition can be relaxed for the nonrelativistic (-Rq ^ 1) long-wavelength {k <^ 1) 
limit to the following instability condition 



e < er 



H 



\ 



Ro 
3 



/3 



Ro 



'1 + Rt, 



(14) 



Thus, for every plasma composition (every value of Z) there exists a -Ro-region in Ro-e plane 
for a given Zeeman parameter e < e^, where the instability sets-in for the magnetosonic 
linear waves and for e > the linear waves are completely stable. This effect has been 
discussed in detail in the next section, where, it has been shown that the width of instability 
region only loosely depends on the wavenumber, k (e.g. see Fig. 4). 



IV. NUMERICAL RESULTS AN ANALYSIS 



It is clearly observed from Eq. (ITT]) that, in the nonrelativistic degeneracy limit of our 
model the diffraction parameter, H, is negligibly small, hence, the magnetosonic waves 
are almost dispersionless, while, for the relativistic degeneracy regime the values of, H 
are relatively large and the linear waves become dispersive. The dispersion in this regime 
becomes more pronounced for higher values of wave number, k. In Fig. 1, the variation 
of paramagnetic susceptibility Xp: given by Eq. @, versus the normahzed plasma mass- 
density, po/Pc is shown. It is observed that the increase in the mass-density (equivalently 
increase in the relativistic degeneracy parameter of the plasma) causes the increase in the 
relativistic susceptibility value. Figure 2 shows the variation of the magnetosonic wave 
dispersion relation due to variations in the relativistic degeneracy parameter, Ro (presenting 
the normalized plasma mass-density) and the normalized Zeeman energy, e (presenting the 
magnetic field-strength) in the two distinct nonrelativistic and relativistic degeneracy limits 
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in the absence of Coulomb interactions {(3 = 0). From Fig. 2(a) it is apparent that in the 
nonrelativistic regime the speed of hnear magnetosonic waves decrease with increase in the 
relativistic degeneracy parameter for a given fixed value of the magnetic field parameter, e. 

On the other hand, it is observed from Fig. 2(b) that (in the absence of Coulomb 
interaction) increase in the strength of the magnetic field leads to the increase in the wave 
speed in this regime. However, as it is revealed from Fig. 2(c), in the relativistic degeneracy 
regime, the increase in the normalized plasma mass-density (i.e. increase in the relativistic 
degeneracy parameter, Rq) results in an increase in the linear magnetosonic speed which is 
opposite to the case of nonrelativistic degeneracy case shown in Fig. 2(a). Furthermore, 
increase in the magnetic field strength in relativistic degeneracy regime in ineffective in 
changing the linear wave-speed for a fixed normalized i?o-parameter which is contrary to the 
case of nonrelativistic degeneracy (e.g. compare Figs. 2(b) and 2(d)). 

Distinguished feature is observed in Fig. 2(e) where a minimum in phase speed exists at 
a given value of relativistic degeneracy, Rom- This minimum {Rom) occurs in nonrelativistic 
degeneracy regime value of which depends strongly on the magnitude of the Zeeman energy 
parameter, e. The occurrence of the minimum in wave speed is related to the competition of 
the positive electron degeneracy pressure and the negative pressure-like spin paramagnetism 
effect. As it is revealed comparing Figs. 2(d) and 2(e) the value of Rom increases with 
increase in the strength of the magnetic field. The corresponding value of i?om — 0.1 
(no ~ 5.9 X 10^^ cm^'^) is obtained for e = 10~^ {Bo ~ lO^T), while, this value increases 
to i?om ^ 0.25 (no ~ 9.2 x lO^^cm"^) for e = 10"^ {Bq ~ lO^T). Figure 3(a), clearly 
confirms the above statement. On the other hand. Fig. 3(b) reveals that the value of 
the minimum phase-speed of the linear waves is nearly independent of the value of the 
wave-number, k. The observed differences in dispersion of linear magnetosonic waves in 
different plasma degeneracy regimes mark distinct behavior of linear magnetosonic wave in 
dense paramagnetic quantum plasmas which is typical of the astrophysical situations and 
can improve our understanding of the linear wave propagation mechanism in superdense 
highly magnetized collapsing white dwarfs or neutron stars. Such distinctive features has 
already been suggested to be present in nonlinear wave dynamics of ion-acoustic solitons in 



Thomas- Fermi plasmas 36l. |37|. 

Figure 4(a) shows and compares the variation of dispersion relations for '^He and ^^Fe 
plasma compositions for a fixed value of Zeeman energy parameter. It is revealed that for 
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the ^^Fe magnetized plasma there is a instabihty gap for a range of relativistic degener- 
acy parameter, Rq for which magnetosonic hnear waves can not propagate. The width of 
the instabihty gap seems not to vary significantly with the wavenumber value, k. Figure 
4(b) reveals that for the ^^Fe composition also the instability sets-in when the magnetic 
field strength is lowered. Remarkably it is revealed, comparing Figs. 4(a) and 4(b) that 
the change in the magnetic field strength for a fixed plasma composition only affects the 
nonrelativistic degeneracy part of the plasma, while, the composition change in plasma for 
fixed magnetic field strength impacts the relativistic degeneracy part. The instability region 
for compositions Z = 2,4, 6, 26 are shown in Fig. 5. It is also remarked that the instability 
area increases with increase in the atomic-number of plasma composition. As it is evident 
from Fig. 5, for every composition there is a minimum Zeeman parameter, e, below which 
the instability of linear magnetosonic waves occur. For a given e- value below this minimum 
the instability is bounded by two values of degeneracy parameter Rq. The Rq range for 
instability increases as the magnetic field strength is lowered, however, for iron plasma com- 
position the situation is different and the plasma becomes linearly unstable for all magnetic 
field strength values above a minimum i?om value. These features might be crucial on the 
stability of compact stars such as white dwarfs and neutron star crusts. For instance, the 
crust of a neutron star or a white-dwarf core is believed to be iron-rich while they are also 
relativistically degenerate and might be highly magnetized. Also, the outer regions of a 
white-dwarf star may contain degenerate carbon or helium. There has been many compact 
star reports with strong internal or external magnetic fields 381-141 1. 
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FIGURE CAPTIONS 



Figure- 1 

The paramagnetic susceptibility, Xp, versus normalized electron number-density, Uc/Nq, 
for a Fermi- Dirac plasma. 

Figure-2 

(Color online) Variation of the magnetoacoustic dispersion curves with respect to the 
changes in the relativistic degeneracy parameter, Rq and the normalized Zeeman-energy, e, 
in the nonrelativistic plasma degeneracy regime (Figs. 1(a) and 1(b)) and in the relativistic 
plasma degeneracy regime (Figs. 1(c) and 1(d)). The minimum dispersion slope (wave 
speed) is shown for two different values of magnetic field strength (Figs. 1(e), 1(f)). 

Figure-3 

(Color online) Variation of the magnetoacoustic wave phase speed with respect to the 
relativistic degeneracy parameter, Rq for differnt values of the normalized Zeeman-energy, 
e (a), the wave-number, k (b). The dash size of curve increases appropriately as the varied 
parameter in each plot is increased. 

Figure-4 

(Color online) Figure 4(a) shows the complete stabihty for ife^ (top surface) but partial 
stability for ^^Fe (bottom surfaces) plasma composition at e = 10"^. Figure 4(b) shows the 
instability inset for ^^Fe as the Zeeman parameter (magnetic-field strength) decreases from 
e = 10-6 to e = 10-^ 

Figure-5 

(Color online) Unstable linear magnetosonic wave regions in Ro-e plane for different 
plasma compositions, Z — 2,4,6,26 with condition Rq « 1. The dashed areas (respec- 
tively from small to large correspond to Z — 2, 4, 6, 26) indicate the unstable regions for 
different plasma compositions. 



13 



